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Abstract 


We  gathered  and  processed  infrared  (IR)  images  taken  by  a  quantum  well 
infrared  photodetector  (QWIP)  camera  to  support  the  development  of  an 
algorithm  for  airborne  target  detection.  This  report  addresses  the 
statistical  aspects  of  this  effort.  Digitized  images  consisted  of 
background  clutter  alone  and  background  clutter  with  a  small  aircraft 
target.  We  used  a  template-matching  technique  to  detect  the  small 
aircraft  amidst  the  background  clutter.  We  propose  a  reasonable  statistic 
for  use  in  our  algorithm.  Our  primary  focus  is  to  describe  the 
experimental  design  and  analysis  used  in  determining  parameter  settings 
for  the  detection  algorithm  and  to  report  those  results. 
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1.  Introduction 


We  gathered  and  processed  infrared  (IR)  images  taken  by  a  quantum  well 
infrared  photodetector  (QWIP)  camera  to  support  the  development  of  an 
algorithm  for  airborne  target  detection.  This  report  addresses  the  statistical 
aspects  of  this  effort.  Digitized  images  consisted  of  background  clutter  alone 
and  background  clutter  with  a  small  aircraft  target.  We  used  a  template¬ 
matching  technique  to  detect  the  small  aircraft  amidst  the  background  clut¬ 
ter.  We  propose  a  reasonable  statistic  for  use  in  our  algorithm.  Our  primary 
focus  is  to  describe  the  experimental  design  and  analysis  used  in  determin¬ 
ing  parameter  settings  for  the  detection  algorithm  and  to  report  those  re¬ 
sults. 

Multiple  tapes  of  IR  imagery  were  collected.  Two  aircraft  flew  over  the  same 
area  at  different  altitudes.  A  Cessna  aircraft  was  used  as  a  surrogate  target 
vehicle.  The  Cessna's  path  crossed  that  of  the  observing  aircraft,  flying  at  a 
greater  altitude  with  its  IR  sensor  pointing  straight  down.  Each  time  the 
two  aircraft  crossed,  the  Cessna's  flight  was  recorded  by  the  continuous 
QWIP  imagery.  Figures  1  through  5  show  representative  imagery.  One  hun¬ 
dred  fifty-two  images  were  selected  from  these  tapes,  showing  the  Cessna 
at  different  apparent  ranges  and  aspect  angles  and  against  a  variety  of  natural 
background  clutter.  Fifty  additional  images  without  the  Cessna  were  also 
selected.  These  additional  images  included  more  confusing  man-made  clut¬ 
ter.  All  images  were  digitized  with  512  x  512  8-bit  pixels  with  grayscale 
values  ranging  from  white  (255)  to  black  (0). 

Our  goal  was  to  develop  a  reliable  algorithm  for  automatic  detection  of  the 
Cessna. 


Figure  1.  QWIP  IR. 
Image  of  Cessna 
against  natural 
background.  Provided 
by  Naval  Research 
Laboratory. 
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Figure  2.  QWIP IR. 
Image  of  Cessna 
against  natural 
background.  Provided 
by  Naval  Research 
Laboratory 


Figure  3.  QWIP  IR. 
Image  of  Cessna 
against  natural 
background.  Provided 
by  Naval  Research 
Laboratory. 
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Figure  4.  QWIP IR. 
Image  of  man-made 
background.  Provided 
by  Naval  Research 
Laboratory. 


Figure  5.  QWIP  IR. 
Image  of  man-made 
background. 
Provided  by  Naval 
Research  Laboratory. 


2.  Template  Matching 


Before  introducing  the  complete  algorithm,  we  briefly  discuss  template 
matching,  a  key  component  of  that  algorithm.  As  the  name  suggests,  tem¬ 
plate  matching  involves  overlaying  a  template  on  the  image  and  testing  to 
determine  if  the  image  and  the  template  match.  A  match  is  said  to  occur  if 
there  is  sufficient  agreement  between  the  edges  of  a  target  template  and  the 
edges  of  the  image.  Since  the  locations  of  the  template  edges  are  known,  the 
interior  and  exterior  of  the  template  are  completely  defined.  In  the  images 
examined,  the  Cessna  was  invariably  quite  cold.  If  the  template  and  target 
were  exactly  aligned,  we  would  expect  pixel-grayscale  values  on  the  inte¬ 
rior  of  the  template  to  be  small  (darker)  and  relatively  homogeneous.  Pixel- 
grayscale  values  exterior  to  the  template  (in  the  background)  would  tend  to 
be  larger  (lighter)  and  less  homogeneous.  Consider  an  edge  of  the  aligned 
template  and  target.  The  difference  across  this  edge  between  an  exterior 
pixel  value  and  a  neighboring  interior  pixel  value  would  tend  to  greatly 
exceed  zero.  This  property  of  an  edge  serves  as  the  basis  for  our  decision 
regarding  whether  the  image  and  template  match. 

We  wish  to  parlay  this  property  into  a  reasonable  algorithm  to  detect  the 
presence  of  the  Cessna  in  our  image.  A  formal  hypothesis  test  is  not  devel¬ 
oped.  Rather,  we  offer  some  justification  for  one  of  the  test  statistics  tried 
and  then  concentrate  on  determining  optimum  parameter  values  for  our 
algorithm  based  on  a  learning  sample  of  images. 

The  statistic  arises  as  follows.  Thirty-eight  ordered  pairs  are  formed  in  fixed 
relation  to  the  template,  with  each  pair  consisting  of  an  exterior  point  and  a 
neighboring  interior  point.  (Pairs  are  numbered  consecutively  around  the 
template.)  Let  Ey  and  ly  denote  the  grayscale  values  of  the  exterior  and 
interior  points,  respectively.  Ignoring  any  dependence,  consider  the  pixels 
in  two  adjacent  ordered  pairs  as  samples  of  size  2  from  the  interior  and 
exterior  regions.  Equation  (1)  can  be  thought  of  loosely  in  terms  of  rationale 
for  analysis  of  variance  in  that  we  are  attempting  to  compare  variation  across 
groups  with  variation  within  groups.  We  take  several  liberties.  The  last  term 
represents  within  variability,  based  only  on  the  pixels  internal  to  the  tem¬ 
plate.  If  we  substitute  an  absolute  value  operator  for  the  traditional  square, 
then  I  Ij  —  Jy^j  I  is  the  sum  of  absolute  deviations  from  the  sample  mean  (Jy  + 
^/+i)/2.  Doubling  this  value  accounts  for  the  contribution  to  within  varia¬ 
tion  made  by  the  external  pixels  imder  the  assumption  that  the  variation  is 
similar.  For  between  variation,  consider  that  (Ey  +  Ey+i)/2  is  an  estimate  for 
the  external  mean  and  that  (Ey  +  Ey^.i  +  ly  +  ly+j)/ 4  is  an  estimate  for  the  grand 
mean.  The  sum  of  squares  between  is  analogous  to  four  times  the  absolute 
difference  between  these  two  terms,  leaving  I  (Ey  +  Ey^j)  -  (7y  +  ly^j)  I .  Finally, 
because  of  the  interest  in  a  one-sided  alternative  with  external  values  being 
expected  to  exceed  internal  values,  the  absolute  value  is  dropped.  We  are 
left  with  a  heuristically  justified  difference  in  eq  (1),  which  will  assume  large 
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values  when  the  two  pairs  each  straddle  a  template  edge  that  is  aligned 
with  the  target. 

Two  related  statistics  that  we  considered  are  given  as  equations  (2)  and  (3). 
Equation  (2)  is  reasoned  along  the  lines  of  (1),  but  within  variation  is  com¬ 
puted  separately  interior  and  exterior  to  the  template.  Equation  (3)  makes 
only  one  comparison  across  a  potential  edge  and  uses  the  neighboring  y+1®* 
internal  pixel  for  a  within-variation  estimate.  The  multiplication  by  two 
adjusts  the  scale  of  equation  (3)  to  that  of  equations  (1)  and  (2). 

(Ej  - 1,)  +  (E;,i  -  -(IE,.-  E,.,i  I  +  1 1,.  - 1,,,  1 )  (2) 

2(E,.-I,-IJ,-f,,J)  (3) 
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3.  Proposed  Algorithm 


The  previous  discussion  suggests  why  proposed  statistics  might  reveal  an 
edge.  To  complete  the  test  of  template  and  target  alignment,  however,  we 
must  look  at  the  entire  template  boundary  as  a  function  of  scale  and  orien¬ 
tation.  For  fixed  scale  and  orientation,  the  template  boundary  evaluation  is 
supported  by  {[(Ey,  7y),  (Ey+„  Jy+i)]:  j  =  1, 38),  with  the  first  point  reused  for/  = 
39.  Thirty-eight  individual  tests  based  on  (1)  to  (3)  record  a  response  if  the 
statistic  exceeds  some  threshold  level  x.  For  the  edges  of  the  template  to  be 
said  to  match  the  target,  the  collective  responses  for  edges  around  the  tar¬ 
get  must  exceed  a  second  threshold,  T.  These  two  parameters,  X  and  T,  de¬ 
termine  whether  a  proposed  template  aligns  with  the  target,  that  is,  whether 
the  target  is  detected. 

Equally  important  is  determining  what  templates  should  be  proposed.  An 
image  of  512  X  512  pixels  is  searched  for  a  subset  that  collectively  yields 
edges  consistent  with  the  target.  The  subset  is  a  collection  of  38  pixel  pairs 
with  fixed  location  relative  to  a  prospective  center  point  for  the  target.  For 
the  template  to  match  the  target,  both  must  have  approximately  the  same 
angular  orientation  in  the  view  considered.  Thus,  many  orientations  must 
be  tried  to  achieve  coverage  of  possible  orientations.  Additionally,  the  tem¬ 
plate  must  be  adjusted  for  various  ranges.  A  target  overlain  by  a  template 
with  the  true  center  point  and  the  correct  orientation  but  very  different  scale 
will  likely  not  be  detected.  Thus,  a  sufficiently  fine  resolution  of  scales  must 
also  be  taken  into  account.  Considering,  for  example,  each  pixel  as  a  poten¬ 
tial  center  point,  36  orientations  and  8  scales  yields  a  set  of  75,497,472  pos¬ 
sible  templates  to  be  proposed. 

To  reduce  the  number  of  computations,  we  perform  the  following  prepro¬ 
cessing.  At  what  we  will  refer  to  as  level  1,  we  consider  points  as  potential 
template  centers  only  if  their  grayscale  value  <  5.  This  is  reasonable  since 
the  target  center  shows  up  nearly  black  on  the  images  we  have.  For  each  of 
these  center  points,  we  take  a  sparse  sample  of  orientation  x  range.  Azi¬ 
muths  are  taken  every  30°  instead  of  every  10°.  Two  rather  than  eight  scales 
are  tested.  Only  seven  ordered  pairs  about  the  template  are  considered  for 
an  edge  test.  For  each  center  point  passing  this  reduced  edge  test,  we  also 
flag  the  eight  adjacent  neighbors.  (Some  multiple  counting  would  be  ex¬ 
pected.)  Thus,  the  preprocessing  would  result  in  a  set  of  potential  center 
points,  along  with  their  nearest  neighbors,  that  had  passed  the  grayscale 
screen  and  reduced  edge  test  screen  for  some  orientation  and  scale. 

The  proposed  algorithm  detects  a  target  if  a  pixel  passing  through  prepro¬ 
cessing  for  a  specific  orientation  and  scale  results  in  more  than  T  edge  tests 
about  the  template  exceeding  the  threshold  x. 
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4.  Strategy  for  Optimization 


A  response  surface  approach  was  used  to  fix  the  four  parameters  of  the 
model  (x,  T,  number  of  scales,  number  of  orientations)  and  to  determine 
which  of  the  three  test  statistics  performed  best.  Two  responses  are  consid¬ 
ered:  the  maximum  probability  of  detection  while  maintaining  zero  false 
alarms  (Max_P(D))  and  the  area  under  a  receiver  operator  curve  (ROC), 
ROC_Area.  For  each  setting  of  the  algorithm  parameters,  we  evaluated  two 
sets  of  images — ^those  containing  a  target  and  those  not  containing  a  target. 
Consider  an  orientation,  scale,  and  edge  test  threshold  x.  The  ROC  arises  by 
recording  the  proportion  of  images  whose  best  pixel  passes  at  least  0,  1, 
2, ...  ,38  edge  tests.  All  images,  with  or  without  targets,  will  pass  zero  edge 
tests  and  earn  a  proportion  1.  See  figure  6.  Few  images  with  targets  and, 
ideally,  no  images  without  targets  would  pass  all  38  edge  tests  to  earn  pro¬ 
portions  of  approximately  zero.  The  goal  in  examining  the  ROC  is  to  deter¬ 
mine  the  number  of  edge  tests  (0  to  38)  that  make  the  algorithm  sufficiently 
sensitive  to  detect  a  target  without  many  false  alarms.  An  area  vmder  the 
curve  of  near  1  is  indicative  of  such  a  test. 

In  our  computer  experimentation,  we  examine  the  impact  of  orientation, 
scale,  and  X  for  each  of  three  proposed  statistics  using  the  ROC_Area  and 
Max_P(D),  with  an  eye  toward  optimizing  those  responses.  The  threshold 
T  for  the  number  of  edge  tests  required  to  conclude  detection  is  determined 
from  examining  near-optimum  ROCs. 


Figure  6.  A  promising 
ROC. 


S? 

c3 


1.1 

.9 
,7 

.5 

.3 
.1 

-.2  0  .2  A  .6  .8  1  1.2 


Trend  plot 


i 

^  ' 

i\ 

i 

i 

38 

Detect  target  in  a  nontarget  image  (%) 


7 


5.  Learning  Sample  Experiments 


Twenty  images  containing  targets  were  selected.  Portions  of  these  images 
with  no  target  present  were  used  as  nontarget  images  (whole  image  less  a 
60  X  60-pixel  area  about  the  target).  Ten  additional  pure  clutter  images  were 
also  used.  Thus,  the  target  images  were  considered  to  number  20  and  the 
nontarget  images,  30  for  the  training  testing. 

5.1  Learning  Set  Analysis  Phase  I 

The  computer  experiments  for  the  initial  test  region  were  rim  with  the  fac¬ 
tors  and  levels  given  in  table  1.  The  threshold  x  is  the  amount  a  chosen  edge 
test  statistic  must  exceed  for  one  to  conclude  that  an  edge  is  present.  The 
orientation  is  the  number  of  steps  taken  over  360°  for  template  positioning. 
For  example,  24  steps  correspond  to  15°.  Stepping  is  performed  relative  to  a 
random  start  orientation.  Scale  indicates  the  number  of  uniformly  spaced 
divisions  over  the  interval  of  target  ranges  considered.  The  statistic  refers 
to  the  expressions  introduced  in  equations  1  to  3.  We  ran  a  complete  facto¬ 
rial,  using  four-  and  five-way  interactions  as  error. 

Analysis  of  variance  indicated  that  only  the  first  three  factors  of  table  1  were 
significant.  See  the  main  effects  plot  in  figure  7.  Interestingly,  the  statistic 
used  was  not  significant.  No  two-  or  three-way  interactions  were  found  to 
be  significant.  See  figure  8.  This  suggests  that  the  algorithm  parameters  could 
be  adjusted  independently.  A  regression  model  of  ROC_Area  as  a  function 
of  only  the  three  significant  main  effects  explained  79.2  percent  of  the  varia¬ 
tion  and  satisfied  the  usual  regression  analysis  assumptions.  Figures  9  and 
10  show  the  best  and  worst  ROC  over  the  factor  levels  tested.  Figure  9  cor¬ 
responds  to  resolution  on  orientation  of  24,  resolution  for  scale  of  6,  and  x  = 
75,  whereas  figure  10  is  based  on  values  of  8,  2,  and  135  for  orientation, 
scale,  and  x,  respectively. 

Similar  results  were  seen  with  the  use  of  MaxP(D)  as  the  response.  The  same 
three  factors  are  significant  and  a  well-behaved  regression  model  explains 
75.5  percent  of  the  variation. 

Both  responses  showed  room  for  improvement.  The  analysis  suggested  that 
each  might  be  moved  closer  to  optimum  by  (1)  increasing  the  resolution  of 
scale,  (2)  increasing  the  resolution  of  orientation,  and  (3)  lowering  the  thresh¬ 
old  for  detection  of  an  edge. 


Table  1.  Initial  test 
region  factor  levels. 


Factor 

Levels 

T 

75, 90, 105, 120, 135 

Orientation 

24, 12, 8 

Scale 

2,4,6 

Statistic 

1,2,3 

8 


Orientation 


Figure  10.  Worst  ROC. 
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A  second  group  of  experiments  was  run  based  on  what  we  learned  in  the 
first  phase.  Table  2  lists  the  factors  and  levels  used  in  this  second  set  of  tests. 
The  three  statistics  (1)  to  (3)  were  all  retained  for  this  test. 

The  analysis  proceeded  similarly  to  that  of  phase  I.  Again,  only  the  three 
main  effects  drove  the  model.  No  significant  interactions  were  present.  Fig¬ 
ure  11  shows  the  main  effects.  From  this  graph  we  can  see  that  the  ROC_Area 
range  in  which  the  testing  is  done  is  much  higher  than  in  the  first  phase  of 
testing,  all  means  showing  areas  in  excess  of  0.91.  We  also  see  that  the  influ¬ 
ence  of  X  seems  to  level  off  as  the  ROC_Area  nears  0.95.  Values  of  X  less  than 
or  equal  to  60  yield  comparable  results.  A  best-fitting  regression  model  does 
include  a  square  term  for  x,  so  the  apparent  curvature  is  probably  real.  The 
explained  variation  over  this  more  narrow  range  of  responses  is  85.8  per¬ 
cent. 

Boxplots  incorporating  the  data  from  phases  I  and  II  were  constructed  to 
show  the  distribution  of  response  measures  as  a  function  of  the  three  sig¬ 
nificant  factors.  See  figures  12  to  14.  We  wished  to  determine  the  minimum 
scale  and  orientation  and  maximum  X  that  could  be  used  to  provide  a  near 
optimal  ROC.  Algorithm  speed  was  a  concern.  More  scales  and  orienta¬ 
tions  improve  detection  at  the  cost  of  computational  efficiency.  A  larger  dif¬ 
ferential  across  an  edge  was  required  for  x  so  that  we  could  be  sure  that  the 
apparent  edge  was  not  just  a  difference  in  clutter  values.  In  examining  these 
three  figures,  we  conclude  that  there  should  be  8  resolutions  of  scale,  36 
orientations,  and  a  value  of  about  60  for  x.  A  refinement  for  orientation  was 
made  after  looking  at  specific  plots  of  individual  values.  In  our  judgment, 
24  orientations  would  suffice. 

The  parameter  for  our  algorithm  not  yet  established  is  the  threshold  for  the 
number  of  edge  tests  out  of  the  38  that  must  be  passed  for  a  detection  to  be 
said  to  occur.  Our  approach  to  finding  this  threshold  was  to  examine  the 
ROCs  for  which  the  ROC_Area  was  very  high  and  consistent  with  our  opti¬ 
mal  settings  for  the  other  parameters.  For  each  of  these  we  recorded  the 
number  of  edge  tests  passed  corresponding  to  the  Max_P(D).  The  results 
did  not  vary  greatly,  all  being  34, 35,  or  36.  We  chose  to  set  the  threshold  for 
the  number  of  edge  tests  passed  at  35. 
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Table  2.  Phase  II  Test  Factor 

region  factor  levels.  ^ 

Orientation 

Scale 

Statistic 


Figure  11.  Main  effects 
plot. 
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Figure  12.  Boxplots  of 
ROC  Area  by  scale. 
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Figure  13.  Boxplots  of 
ROC_Area  by 
orientation. 
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Figure  14.  Boxplots  of 
ROC_Area  by  r. 
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6.  Test  Sample  Experiment 


A  final  test  was  run  to  confirm  the  choice  of  parameters  for  our  model:  T  = 
60,  orientation  =  24,  scale  =  8,  and  edge  test  threshold  =  35.  We  chose  to  test 
one  level  to  either  side  of  the  optimum.  The  threshold  values  of  34, 35,  and 
36  were  considered  for  each  of  the  27  experimental  combinations  from  a 
factorial  crossing  of  the  factors  in  table  3. 

One  hundred  thirty-two  images  containing  targets  were  selected  for  the 
test  phase.  Portions  of  these  images  with  no  target  present  were  used  as 
nontarget  images  (whole  image  less  a  60  x  60-pixel  area  about  the  target). 
An  additional  40  pure  clutter  images  were  also  used.  Thus  the  target  im¬ 
ages  were  considered  to  number  132  and  the  nontarget  images  172  for  the 
test  phase.  These  images  were  similar  but  not  identical  to  those  used  in  the 
previous  phases. 

Figure  15  shows  the  results  of  the  27  trials  formed  from  the  factor  crossings 
of  table  3.  Results  are  given  in  terms  of  the  two  responses,  ROC_Area  and 
Max_P(D).  As  expected,  the  chosen  parameter  settings  produced  results  in 
the  midrange  for  each  response  over  the  27  trials  for  the  larger  test  data  set. 
All  of  the  runs  over  the  test  data  set  are  included  in  table  4.  It  is  clear  from 
scanning  this  table  that  all  experimentation  is  in  a  region  showing  good 
results  for  both  response  measures.  The  ideal  tradeoff  between  probability 
of  detection  and  computational  efficiency  will  be  dependent  upon  the  ap¬ 
plication  in  which  it  is  used. 


Table  3.  Test  phase 

Factor 

Levels 

factor  levels. 

T 

50, 60, 70 

Orientation 

36, 24, 18 

Scale 

6, 8, 10 

Figure  15.  Algorithm 
performance. 
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Table  4.  Summary 
algorithm  performance. 


Max_P(D) 

ROC_Area 

Scale 

Orientation 

X 

Threshold 

0.962 

0.984 

6 

36 

50 

36 

0.955 

0.976 

6 

36 

60 

35 

0.939 

0.971 

6 

36 

70 

34 

0.955 

0.984 

6 

24 

50 

36 

0.932 

0.976 

6 

24 

60 

35 

0.917 

0.969 

6 

24 

70 

35 

0.939 

0.981 

6 

18 

50 

36 

0.939 

0.976 

6 

18 

60 

35 

0.917 

0.968 

6 

18 

70 

34 

0.939 

0.984 

8 

36 

50 

37 

0.939 

0.975 

8 

36 

60 

35 

0.932 

0.971 

8 

36 

70 

35 

0.955 

0.983 

8 

24 

50 

36 

0.932 

0.975 

8 

24 

60 

35 

0.932 

0.972 

8 

24 

70 

34 

0.955 

0.984 

8 

18 

50 

36 

0.939 

0.975 

8 

18 

60 

35 

0.902 

0.969 

8 

18 

70 

34 

0.962 

0.984 

10 

36 

50 

36 

0.939 

0.976 

10 

36 

60 

35 

0.924 

0.970 

10 

36 

70 

35 

0.955 

0.983 

10 

24 

50 

36 

0.939 

0.975 

10 

24 

60 

35 

0.917 

0.969 

10 

24 

70 

35 

0.955 

0.979 

10 

18 

50 

36 

0.932 

0.975 

10 

18 

60 

35 

0.924 

0.971 

10 

18 

70 

34 
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7.  Discussion 


In  this  paper,  we  used  a  mix  of  engineering  practice  and  statistics  to  try  to 
arrive  at  an  algorithm  for  airborne  target  detection.  Admittedly,  the  approach 
was  somewhat  informal  from  a  statistical  perspective,  but  it  seemingly  pro¬ 
duced  good  results.  Extension  to  other  clutter  backgrounds  would  require 
additional  testing. 

We  did  consider  more  traditional  approaches,  specifically  a  paired-t  signed 
rank  test  and  a  Wilcoxon  matched-pairs  signed  rank  test.  In  each  instance, 
we  must  overlook  the  spatial  correlation  present  among  neighboring  inte¬ 
rior  and  neighboring  exterior  points.  Threshold  T  is  analogous  to  the  sig¬ 
nificance  level  of  the  test.  Threshold  T  is  analogous  to  a  general  mean  differ¬ 
ence  to  be  exceeded.  Setting  these  values  with  respect  to  our  optimization 
goals  requires  reasoning  along  the  lines  used  in  this  paper.  With  assump¬ 
tions  for  these  procedures  in  question  and  with  no  real  advantage  to  the 
general  modeling  approach,  we  settled  on  the  approach  presented  here. 


8.  Acknowledgment 


The  authors  would  like  to  thank  the  Ballistic  Missile  Defense  Organization 
and  the  Naval  Research  Laboratory  for  providing  the  data  and  this  oppor¬ 
tunity  to  explore  target-detection  algorithms. 


15 


9.  References 


Baras,  J.S.,  and  D.C.  MacEnany  (1992),  "Model  based  ATR  algorithms  based  on 
reduced  target  models,  learning  and  probing,"  Proc.  2"'*  of  the  Second  Au¬ 
tomatic  Target  Recognizer  Systems  and  Technology  Conference,  1992. 

Bienenstock,  D.,  D.  Geman,  S.  Geman,  and  D.E.  McClure  (1990),  Phase  II  Technical 
Report,  Development  of  Laser  Radar  ATR  Algorithms,  Contract  No.  DAAL02- 
89-C-0081,  Night  Vision  and  Electronic  Sensors  Directorate  (CECOM),  Ft. 
Belvoir,  VA. 

Nguyen,  D.M.  (1990),  An  Iterative  Technique  for  Target  Detection  and  Segmenta¬ 
tion  in  IR  Imaging  Systems,  Ft.  Belvoir,  VA. 

Severson,  W.E.  (1996),  The  SSV  Vehicle's  FLIR  Target  Detection  Capability, 
Lockheed  Martin  Astronautics,  Contract  No.  DASG60-95-C-0062,  U.S.  Army 
Space  &  Strategic  Defense  Command. 


16 


Distribution 


Adirmstr 

Defns  Techl  Info  Ctr 

ATTN  DTIC-OCP 

8725  John  J  Kingman  Rd  Ste  0944 

FT  Belvoir  VA  22060-6218 

DARPA 

ATTN  S  Welby 
3701  N  Fairfax  Dr 
Arlington  VA  22203-1714 

Ofc  of  the  Secy  of  Defns 
ATTN  ODDRE  (R&AT) 

The  Pentagon 

Washington  DC  20301-3080 

AMCOM  MRDEC 

ATTN  AMSMI-RD  W  C  McCorkle 

Redstone  Arsenal  AL  35898-5240 

US  Army  TRADOC 

Battle  Lab  Integration  &  Techl  Dirctrt 

ATTN  ATCD-B 

FT  Monroe  VA  23651-5850 

BMDO 

ATTN  L  Lome 

7100  Defense  Pentagon 

Washington  DC  20301-7100 

Hdqtrs  Dept  of  the  Army 
ATTN  DAMO-FDT 
400  Army  Pentagon  Rm  3C514 
Washington  DC  20301-0460 

US  Military  Acdmy 
Mathematical  Sci  Ctr  of  Excellence 
ATTN  MADN-MATH  MAJM  Johnson 
Thayer  Hall 

West  Point  NY  10996-1786 
Dir  for  MANPRINT 

Ofc  of  the  Deputy  Chief  of  Staff  for  Prsnnl 
ATTN  J  Hiller 
The  Pentagon  Rm  2C733 
Washington  DC  20301-0300 

SMC/CZA 

2435  Vela  Way  Ste  1613 
El  Segundo  CA  90245-5500 


US  Army  ARDEC 
ATTN  AMSTA-AR-TD 
Bldgl 

Picatinny  Arsenal  NJ  07806-5000 

US  Army  Info  Sys  Engrg  Cmnd 
ATTN  AMSEL-IE-TD  F  Jenia 
FT  Huachuca  AZ  85613-5300 

US  Army  Natick  RDEC  Acting  Techl  Dir 
ATTN  SBCN-T  P  Brandler 
Natick  MA  01760-5002 
US  Army  Simulation  Train  &  Instrmntn 
Cmnd 

ATTN  AMSTI-CG  M  Macedonia 

ATTN  J  Stahl 

12350  Research  Parkway 

Orlando  FL  32826-3726 

US  Army  Tank-Automtv  Cmnd  RDEC 

ATTN  AMSTA-TR  J  Chapin 

Warren  MI  48397-5000 

Nav  Rsrch  Lab 

ATTN  Code  7604  J  Goldspeil 

4555  Overlook  Ave  SW  Bldg  209  Rm  124 

Washington  DC  20375-5352 

Nav  Rsrch  lab  Sensor  &  Data  Analysis  Sect 
ATTN  Code  5621  F  Bucholtz 
Washington  DC  20375-5000 

Nav  Surfc  Warfare  Ctr 
ATTN  CodeB07  J  Pennella 
17320  Dahlgren  Rd  Bldg  1470  Rm  1101 
Dahlgren  VA  22448-5100 

Univ  of  Texas  at  Austin 
Inst  for  Advncd  Tchnlgy 
PO  Box  202797 
Austin  TX  78720-2797 

Hicks  &  Assoc  Inc 
ATTN  G  Singley  III 
1710  Goodrich  Dr  Ste  1300 
McLean  VA  22102 


17 


Distribution 


US  Army  Rsrch  Lab 
ATTN  AMSRL-CI-CT  A  Neiderer 
ATTN  AMSRL-CI-CT  B  Bodt  (5  copies) 
ATTN  AMSRL-CI-CT  C  Hansen 
ATTN  AMSRL-CI-CT  E  Heilman 
ATTN  AMSRL-CI-CT  G  Moss 
ATTN  AMSRL-CI-CT  J  Forester 
ATTN  AMSRL-CI-CT  J  O'May 
ATTN  AMSRL-CI-CT  M  Thomas 
ATTN  AMSRL-CI-CT  P  Jones 
ATTN  AMSRL-CI-CT  RC  Kaste 
ATTN  AMSRL-CI-LP  (305) 

Aberdeen  Proving  Ground  MD  21005 


Director 

US  Army  Rsrch  Lab 

ATTN  AMSRL-RO-D  JCI  Chang 

PO  Box  12211 

Research  Triangle  Park  NC  27709 
US  Army  Rsrch  Lab 

ATTN  AMSRL-CI-CB  P  David  (3  copies) 
ATTN  AMSRL-CI-CN  D  Hillis  (5  copies) 
ATTN  AMSRL-CI-IS-R  Mail  &  Records  Mgmt 
ATTN  AMSRL-CI-IS-T  Techl  Pub 
ATTN  AMSRL-CI-OK-TL  Techl  Lib  (2  copies) 
ATTN  AMSRL-D  D  R  Smith 


18 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1 21 5  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Papenvork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1 .  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

August  2001  Final  FY98 

4.  TITLE  AND  SUBTITLE  Att  AppHcation  of  Statistics  to  Algorithm 

Development  in  Image  Analysis 

5.  FUNDING  NUMBERS 

DA  PR:  AH48 

PE:  611102H48 

6.AUTHOR(S)  Barry  Bodt,  Philip  David,  and  David  Hillis 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 

Attn:  AMSRL-CI-CD  email:  babodt@arl.army.mil 

2800  Powder  Mill  Road 

Adelphi,  MD  20783-1197 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

ARL-TR-2333 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Ballistic  Missile  Defense  Organization 

7100  Defense,  Pentagon 

Washington,  DC  20301-7100 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

1 1 .  SUPPLEMENTARY  NOTES 

ARL  PR:  IFEPCl 

AMS  code:  611102H4811 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT  Approved  for  pubUc  release;  distribution 
unlimited. 

12b.  DISTRIBUTION  CODE 

1 3.  ABSTRACT  (Maximum  200  words) 

We  gathered  and  processed  infrared  (IR)  images  taken  by  a  quantum  well  infrared  photodetector 
(QWIP)  camera  to.  support  the  development  of  an  algorithm  for  airborne  target  detection.  This 
report  addresses  the  statistical  aspects  of  this  effort.  Digitized  images  consisted  of  background 
clutter  alone  and  background  clutter  with  a  small  aircraft  target.  We  used  a  template-matching 
technique  to  detect  the  small  aircraft  amidst  the  background  clutter.  We  propose  a  reasonable 
statistic  for  use  in  our  algorithm.  Our  primary  focus  is  to  describe  the  experimental  design  and 
analysis  used  in  determining  parameter  settings  for  the  detection  algorithm  and  to  report  those 
results. 

14.  SUBJECT  TERMS  infrared  sensor,  target  detection,  experimental  design 

IS.  NUMBER  OF  PAGES 

25 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION 

OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified  Unclassified 

20.  LIMITATION  OF  ABSTRACT 

UL 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39-18 


298-102 

19 


